Can a Bose gas be saturated? 
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Bose-Einstein condensation is unique among phase transitions between different states of matter in the sense 
that it occurs even in the absence of interactions between particles. In Einstein's textbook picture of an ideal 
gas, purely statistical arguments set an upper bound on the number of particles occupying the excited states of 
the system, and condensation is driven by this saturation of the quantum vapour |1|. Dilute ultracold atomic 
gases are celebrated as a realisation of Bose-Einstein condensation in close to its purely statistical form |2}{6l. 
Here we scrutinise this point of view using an ultracold gas of potassium (^^K) atoms, in which the strength of 
interactions can be tuned via a Feshbach scattering resonance III (8). We first show that under typical experi- 
mental conditions a partially condensed atomic gas strongly deviates from the textbook concept of a saturated 
vapour. We then use measurements at a range of interaction strengths and temperatures to extrapolate to the 
non-interacting limit, and prove that in this limit the behaviour of a Bose gas is consistent with the saturation 
picture. Finally, we provide evidence for the universality of our observations through additional measurements 
with a different atomic species, ^^Rb. Our results suggest a new way of characterising condensation phenomena 
in different physical systems. 



Einstein's prediction of a purely statistical phase transition 
follows directly from the Bose distribution for an ideal (non- 
interacting) gas of identical particles |[5l|6l[9l. At a given tem- 
perature T, there is an upper bound A^c^^^ (T) on the number 
N' of bosons in the excited states of the system. For the ex- 
perimentally relevant case of a 3D harmonically trapped gas 
such that k^T/Ti is large compared to the trapping frequen- 
cies uji along each axis {i = z), the saturation of the gas 
is given by ifTOl 

TV' < iV('d) = C(3) , (1) 

where Co is the geometric mean of the cj^'s and ( is the Rie- 
mann function [C(3) ~ 1.202]. Increasing the total num- 
ber of particles A/tot above the critical value A'c^^^ results in 
the macroscopic occupation of the ground state, i.e. Bose- 
Einstein condensation. 

The condensation observed in weakly interacting atomic 
Bose gases is generally believed to provide a faithful illus- 
tration of the statistical phase transition proposed by Einstein. 
At the same time, differences from ideal gas condensation are 
also observed lfTTUT4l . for example in the small deviations of 
the measured critical atom number Nc from the ideal gas pre- 
diction A^i^"^^ CBII3. Here we focus on the concept of satura- 
tion as the underlying mechanism driving the transition. One 
might expect that the saturation inequality ([T]) is essentially 
satisfied in these systems, with just the value of the bound on 
the right-hand side slightly modified. We prove that this is 
far from being the case, and show how to reconcile experi- 
mental findings with the prediction ([T]). The crucial step in 
our work is a proper disentanglement of the subtle role of in- 
teractions in condensation. While it is correct that Einstein's 
statistical argument does not explicitly invoke interactions be- 
tween the particles, it does assume that the gas is in thermal 
equilibrium, which is fundamentally impossible to attain in 



a completely non-interacting system ifTS l. We overcome this 
problem by making measurements at a range of interaction 
strengths, always sufficient to ensure thermal equilibrium, and 
then extrapolating our results to the non-interacting limit. 

We perform conceptually simple experiments in which we 
keep the temperature of the gas constant and vary the atom 
number. Our experimental scheme is outlined in Fig. [T^ (see 
also Supplementary Information). We start with a partially 
condensed gas of ^^K atoms in the |F, m^) = |1, 1) lower 
hyperfine ground state, produced in a crossed optical dipole 
trap |17|. The optical potential near the bottom of the trap 
is close to harmonic, with Co /27r varying between 60 Hz and 
80 Hz for data taken at different temperatures. We tune the 
strength of repulsive interactions in the gas, characterised by 
the positive s-wave scattering length a, by applying a uniform 
external magnetic field in the vicinity of a Feshbach scattering 
resonance centred at 402.5 gauss 1 18 1. We always prepare the 
condensed gas at a = 135 ao, where ao is the Bohr radius, 
and then adjust the scattering length to the desired value by 
changing the applied magnetic field. In a given experimental 
series, the temperature is kept constant by fixing the depth of 
our optical trap, and the atom number is varied by holding the 
gas in the trap for a variable time up to several tens of sec- 
onds. During this time the total atom number slowly decays 
due to three-body recombination, scattering of photons from 
the trapping laser beams and collisions with the background 
gas in the vacuum chamber, while elastic collisions among 
the trapped atoms ensure equilibrium redistribution of parti- 
cles between the condensate and the thermal gas. 

An example of an experimental series, taken at a = 135 ao 
and T = 177 nK, is shown in Fig. [T]3. For each hold time 
between 1 and 110 seconds we extract the number of atoms 
in the condensate, A^o, and in the thermal gas, N\ from a bi- 
modal fit to the density distribution of the gas after 1 8 ms of 
free time-of-flight expansion from the optical trap. We plot 
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FIG. 1 : Lack of saturation of a quantum degenerate atomic Bose 
gas. (a) Schematic of the experimental setup. Potassium (^^K) atoms 
are confined in an optical dipole trap formed by intersecting two in- 
frared laser beams. The temperature of the gas is fixed by the depth 
of the trapping potential, and the strength of interactions is controlled 
by a uniform magnetic field produced by a pair of coils in Helmholtz 
configuration. The number of thermal (N') and condensed (A^o) 
atoms is extracted from the density distribution after free expansion 
from the trap (see Methods), (b) N' (red points) and A^o (blue points) 
versus the total atom number A^tot at a temperature T = 177 nK and 
a scattering length a — 135 ao- The corresponding theoretical pre- 
dictions for a saturated gas are shown by red and blue solid lines. 
The critical point A^tot = A'c is marked by a vertical dashed line. 

A^o and N' versus the total number of atoms, A^tot, which is 
extracted independently by a direct summation over the den- 
sity distribution. We find that the relationship A^tot = A^o+^' 
is satisfied for all data points to within 0.5%. The standard 
deviation of the temperature for all the points where the con- 
densate is present is 3 nK (see Supplementary Information). 

The predictions for the number of condensed and thermal 
atoms in a saturated gas are shown in Fig. [TJ) by the blue and 
red solid lines, respectively. Specifically, for A^tot > the 
thermal atom number N' remains constant and equal to A^c- 
The deviation of the experimental data from this prediction is 
striking. As the total number of atoms is increased from the 
measured critical value, ^ 200,000, to 450,000, only half 
of the additional atoms accumulate in the condensate. 

In Fig.|2]we show the results of 18 experimental series taken 
at a wide range of scattering lengths, 40 ao < a < 356 ao, and 
temperatures, 115nK < T < 284 nK. Here we focus on the 
regime A^tot > where the condensate is present, and plot 
Nq versus A^tot — ^c- The solid line shows the prediction for 
a fully saturated thermal gas, A^o = ^tot — ^c- The deviation 
of the data from this prediction is clearly observable in all the 
series, and grows with both a and T. 

To explore the relationship between the non- saturation of 
our Bose gases and the interatomic interactions, we start by 
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FIG. 2: Deviation from the saturation picture at a range of inter- 
action strengths and temperatures. We plot A'o versus Ntot — A'c 
for 18 experimental series, each at fixed a and T. The values of the 
scattering length (40-356 ao) and the temperature (1 15-284 nK) are 
encoded in the colour of the data points. The solid line is the predic- 
tion for a saturated gas, A'o = A^tot — A'c- 

identifying the relevant interaction energy. Due to the large 
ratio between the average densities of the condensed and ther- 
mal fractions, the non-ideal behaviour of the thermal compo- 
nent primarily results from its interaction with the condensate. 
The relevant energy scale is then provided by |[T9ll : 



2 V aho 
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where aho = {Ti/muoY^'^ is the spatial extension of the ground 
state of the harmonic oscillator of frequency cD, and m is the 
atomic mass. The energy /io is the mean-field prediction for 
the chemical potential of a gas with Nq atoms at zero temper- 
ature and in the Thomas-Fermi limit 1 19 |. One can understand 
its scaling with A^o in the following way: For a harmonically 
trapped condensate the potential energy per particle scales as 
(X R^, where R is the characteristic size of the cloud, and 
the mean-field interaction energy per particle is proportional 
to the density, Ei oc Nq/R^. Minimising the sum of these 

1/5 2/5 

two energies we get R oc Nq^ and /io oc • 

Guided by this scaling, in Fig.|3]we plot the thermal atom 
number as a function of A/'q^^, for the same experimental 
series as shown in Fig.[TJ). From here we proceed in two steps: 
First, we show that the initial linear increase ofN' withA^o/^ 
can be quantitatively accounted for by the mean-field Hartree- 
Fock (HF) theory for a harmonically trapped gas. Second, 
for the regime of larger condensates, where the theory does 
not fully reproduce the experimental data, we adopt a more 
heuristic approach that still allows us to prove the concept of 
a saturated gas in the non-interacting limit. 

In the HF approach one treats the thermal fraction as an 
ideal gas, but takes into account repulsive interactions with the 

r(id) 



condensate. Within this theory (1911201 one gets Nc = 
and can predict a linear variation of N' /Nc with the small 



3 



320 



280- 



240 



200 




FIG. 3: Quantifying the lack of saturation. Here N' is plotted as 
a function of A^q ^ for the same series as in Fig. II d. The horizontal 
dotted line is the saturation prediction N' = Nc. The blue line is the 
mean-field Hartree-Fock result for a harmonically trapped gas (see 
text), with a slope Suf = 699. The red line is a linear fit to the data 
in the range corresponding to 0.1 < /io/feT < 0.3, which gives a 
non- saturation slope 5' = 1283db84. The solid black line is a guide 
to the eye based on a second-order polynomial fit. The initial slope 
of this line is indistinguishable from HF theory. 



parameter jiQ / k^T: 



= 1 
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with a = C(2)/C(3) ^ 1.37 (see Supplementary Informa- 
tion). The origin of this non- saturation can be qualitatively 
understood by noting that interactions with the condensate 
modify the effective potential seen by the thermal atoms from 
a parabola into a "Mexican hat" shape; this allows the ther- 
mal component to occupy a larger volume, which grows with 
increasing Nq [21 J. 

From Eqs. ^ and Q we define the non- saturation slope 
5'hf = dN'/d{Nl^^) = 1.37X, where X is a dimen- 
sionless parameter of the form X = ^T^a^/^, with ^ = 
0.5 C(3) 15^/^ {kB/ Tiuofa'^^^. The blue line in Fig.jijcorre- 
sponds to this prediction, with the intercept fixed by the mea- 
sured A/'c. It agrees with the data very well for small conden- 
sates, with A^o ^ 10^. corresponding to /io / ^ 0.1. 

To quantitatively test the prediction of Eq. ^ we took sev- 
eral series at different scattering lengths (a = 56 — 274 ao) 
and temperatures (T = 177 — 317 nK), specifically focusing 
on very small values of Nq. We turn off interactions during 
time-of-flight, so that the small condensates almost do not ex- 
pand and can be reliably detected and characterised in absorp- 
tion imaging. For each series we fit the initial non- saturation 
slope, ^0 = dN' / d{NQ^^) for A^o 0, and compare the re- 
sult with the prediction iShf = 1.37X. As shown in Fig.|4j 
the experiment and theory agree within a few percent. 

The agreement of experiments with Eq. ^ for small Nq is 
the first main quantitative result of this paper, and allows us to 
deduce that the initial non- saturation slope would indeed 



vanish in the non-interacting limit, where /io ^ for any A^o- 
This however does not complete our experimental proof since 
this theory works very well only for small condensed fractions 
(see Fig. [3]). For the larger, and experimentally more typical 
values of Nq, we observe an even more pronounced lack of 
saturation of the thermal gas. 

To quantitatively study non- saturation effects at larger Nq 
we take the following heuristic approach: Although the ob- 
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served increase of N' with Nq ' is not perfectly linear, over a 
broad experimentally relevant range it can be well quantified 
by a coarse-grained slope S = A[A^']/A[A^q^^], as indicated 
in Fig.[3]by the red solid line. In order to treat all the data taken 
at various values of a and T equally, for each experimental se- 
ries we consider the same range of values of /io / ^b^, from 
0.1 to 0.3. Note that this range covers more than an order of 
magnitude of A^o values, and encompasses the bulk of the data 
shown in Fig.|2] 

In Fig. [4] we summarise the non- saturation slopes 5'(a, T) 
for the same 18 experimental series shown in Fig.|2] Within 
experimental error all data points fall onto a straight line when 
plotted against the dimensionless X = ^T^a^/^, supporting 
the assumption that we can still use fiQ/k^T as the relevant 
interaction parameter. To further validate our approach, we 
took additional data with a different atomic species, ^^Rb in 
the |F, mi?) = |2,2) state. In this case a = 99 ao is not 
tuneable, and the two experimental series were taken at tem- 
peratures of 175 nK and 203 nK. These two points are in close 
agreement with the ^^K data. 

The prediction ([T]) for an ideal gas (a = 0) corresponds to 
5* = at the origin of the graph in Fig.|4] The mathematical 
limit X ^ is reached in two physically very different limits, 
a ^ and T ^ 0. In the mundane T ^ limit inevitably 
vanishes for any value of Nq, so trivially 5* = 0. It is therefore 
essential that our experiments show that S depends only on the 
parameter X ex T^a^/^, allowing us to deduce its value in the 
a ^ limit for any fixed T. The solid line in Fig.|4]shows the 
result of an unconstrained linear fit to the data, which gives an 
intercept consistent with zero, S{0) = —20 ± 100. Together 
with the success of HF theory in the small A^o regime, this 
confirms the concept of a saturated Bose gas for a broad range 
of experimentally relevant parameters. 

We expect our results to be generic to experiments on har- 
monically trapped 3D Bose gases with (relatively weak) short- 
range s-wave interactions. However a number of important 
questions remain open. Here we studied only the global prop- 
erties of the gas, inferred from the number of atoms in con- 
densed and thermal components; in the future it would be 
of great interest to also study saturation at the level of lo- 
cal densities of the two components, and effectively measure 
the equation of state for a bulk system |23|. Local saturation 
should depend only on the interactions between particles, but 
could for example be different in systems with very strong 
or long-range interactions, where the emergence of a roton 
minimum significantly increases the density of states for low 
energy excitations. Global saturation additionally depends on 
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FIG. 4: Saturation in the non-interacting limit. The non- saturation 
slopes So and S are plotted versus the dimensionless parameter 
X — ^T'^o?^^ (see text). The 5*0 data is directly compared with 
Hartree-Fock theory, 5'hf = 1.37X (blue line), with no free param- 
eters. For the S data, a linear fit (red line) gives dS/ dX = 2.6 ± 0.3 
and an intercept 5'(0) = —20 di 100, consistent with complete satu- 
ration in the ideal gas limit. The S data is based on the eighteen ^^K 
series shown with the same symbol code in Fig. [2] and two additional 
series taken with ^^Rb (black squares). All vertical error bars are sta- 
tistical. The systematic uncertainty in atom numbers N' and A^o is 
< 10%, corresponding to < 6% uncertainty in 5*0 and S values. The 
horizontal error bars include the 3 Hz uncertainty in the trapping fre- 
quencies uj/2n and (for potassium) the 0.1 gauss uncertainty in the 
position of the Feshbach resonance. 

the external potential, and from purely geometric arguments 
we expect the non-saturation effects to grow with the dimen- 
sionality of the system. In this respect it would be particularly 
interesting to study them for atomic gases in disordered poten- 
tials 1241 . where the possible fractal nature of the fluid shape 
can lead to a non-integer effective dimensionality. 

In conclusion, our work shows that the purely statistical 
picture of a saturated Bose gas is not realised in experiments 
with harmonically trapped atomic vapours. However, extrapo- 
lation of our results to the strictly non-interacting limit allows 
us to confirm the textbook picture of Bose-Einstein conden- 
sation as a purely statistical phase transition in an ideal gas. 
Ultracold atomic gases offer great experimental flexibility for 
further studies of (non-)saturation effects, which may provide 
a fruitful way of classifying different geometries and interac- 
tions in many-body systems. Better understanding, and ulti- 
mately control of these effects could also be of practical im- 
portance. Maximising the purity of the condensate for a given 
attainable temperature and atom number is indeed crucial for 
minimising the decoherence processes in atom interferometry 
and quantum information processing. 

Methods Summary 

The bimodal fit function used to determine N' , Nq and T 
is the sum of the Thomas-Fermi profile for the condensate 
and the dilogarithm function g2 for the thermal component, 
with the chemical potential set to zero in the argument of 



92 |[l3][25l. The fit to the whole distribution is improved by: 
(1) refitting A^' by excluding the central region of the image 
occupied by the condensate, (2) refitting A'o with fixed param- 
eters of the g2 function, and (3) fitting the temperature to the 
high energy wings of the distribution by excluding the central 
region of the cloud corresponding to one thermal radius Rt • 
Small effects of varying the exclusion region between 0.8 Rt 
and 1.2 Rt are included in the error bars. 

Within a given series, measurements with different hold 
times are taken in random order. We sometimes observe small 
(< 10%) temperature drifts as a function of A'o or the time 
through the series. We analytically correct for these drifts 
using the temperature scaling given by the Hartree-Fock the- 
ory. This amounts to plotting N'{T/T'f versus nI^^ {T/T'), 
where T' is assigned to each data point and T is the nominal 
temperature for a series, obtained in the A^o ^ limit. 
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FIG. 5: Experimental procedure, (a) Relevant part of the experimental sequence. The orange and green lines show the value of the optical 
dipole trap (ODT) power and Feshbach field respectively as a function of time. The atoms are cooled by an exponential reduction of the ODT 
power. The power is then ramped up slightly over 0.5 s in order to maintain a steady temperature during the hold time. During the cooling 
the Feshbach field is held constant to give a scattering length of 135 ao- Once the final ODT power is reached the Feshbach field is ramped to 
its final value. The range of ODT powers and Feshbach field values used for different experimental series is shown by the dashed lines, (b) 
Total atom number A^tot versus hold time for the series displayed in the main paper Figs. 1 and 3. The data is fitted with a double exponential 
to guide the eye. (c) Temperature versus hold time for the same series. The temperature for A^tot > Nc shows no drift and has a standard 
deviation of 3nK. 

Hartree-Fock Theory 

The simplest theoretical framework to address the effects of atomic interactions in a Bose-condensed gas is the Hartree-Fock 
approximation |[T|. In this model one treats the thermal atoms as a "non-interacting" gas of density n'(r) confined in the self- 
consistent mean-field potential 2^(no(r) + n'(r)), where g = ^uTi^a/m and no(r) is the condensate density. Two further 
approximations take advantage of the small density of the thermal component: (i) One neglects the influence of the thermal 
atoms on the condensate spatial distribution, which is then given by no(r) = Max[(/io — ^(i*))/ 9-> 0] where /io is defined in Eq. 
(2) of the main paper. This is valid as long as A^o ^ <^ho/<^ (i-e. Nq ^ 100 to 1000 for our parameters), (ii) One neglects the 
mean-field energy 2gn' (r) of the thermal atoms, so that the confining potential for these atoms is simply a 3D "Mexican hat" 
potential given by 

Kff(r) = V{v) + 2^no(r) = \V{v) - /io| + Mo • (4) 

The latter approximation amounts to neglecting the shift of the critical temperature due to interactions. 

The number of thermal atoms in the presence of a condensate with A^o atoms can now be estimated using the semi-classical 
version of the Bose-Einstein distribution: 

^ ' Feff(r)-/ioll -l| d'Td'p. (5) 



^' (2^nf J r""^ L'^bT \2m 
The result at first order in ji^/k^T reads (| 1 1, Eqs. (120) and (122)): 



^,3 I C(2),2 MO ... 

where t = T/Tc^^\ and /cbT^^^^ = huj{N/({3)y^^ is the ideal-gas result for the condensation temperature of a Bose gas with 
A^ atoms. Using N/Nc^-' = we can recast the above equation into 

^' -l + g^T^- (7) 



iVi'<i) C(3) kBT 
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